Thermal conductance in a spin-boson model: Cotunneling and low temperature 

properties 
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Bosonic thermal transport through a two-level system is analyzed at temperatures below and com- 
parable to the two-level energy splitting. It is shown that in the low-temperature regime transport 
is dominated by correlated two-boson processes analogous to electron cotunneling in quantum dots 
under Coulomb blockade. We present a detailed analysis of the sequential-cotunneling crossover and 
obtain essentially an analytic description of the transport problem. Perturbative analysis is com- 
plemented by employing scaling properties of the Ohmic spin-boson model, allowing us to extract 
an anomalous low temperature scaling of thermal conductance. 
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I. INTRODUCTION 

The spin-boson model describes the interaction be- 
tween bosonic modes and a two-level quantum system^ 
It has been employed in a wide variety of phenomena 
exhibiting dissipation and loss of quantum coherence. 
Lately the model has found applications in a nonequi- 
librium setting where two bosonic reservoirs with dif- 
ferent temperatures couple through the two-level system 
(TLS). 2 ~ 9 These considerations are motivated by efforts 
to understand thermal properties of artificial molecular 
and nanoscale structures. The spin-boson model is ap- 
plicable to a generic situation where bosonic fluctuations 
couple through a nonlinear region at low enough temper- 
atures so that a truncation to the lowest two energy lev- 
els is sufficient. The spin-boson model can be regarded 
as a bosonic counterpart of the Anderson model which 
describes dynamics of a single spin-degenerate electron 
orbital coupled to leads. Although these model systems 
provide a highly simplified picture of real physical sys- 
tems, they often capture essential properties successfully. 

The spin-boson model is interesting from the point 
of view of heattronics, the manipulation and control of 
heat in artificial nanostructures i 10 ' 11 Nonlinear elements, 
a two-level system being the simplest, with tunable pa- 
rameters are basic building blocks for manipulation of 
heat flow. Possible realizations of the spin-boson thermal 
transport are, for example, molecular junctions coupling 
two phonon reservoirs or electric circuits coupled through 
a superconducting qubit ^2 see Fig. [TJ In the latter ex- 
ample the qubit can be manipulated by external fields 
allowing a flexible control of its properties and enabling 
the investigation of different transport regimes. 

In this paper we present a perturbative theory of ther- 
mal transport in the spin-boson model for arbitrary bath 
temperatures. At temperatures above or comparable to 
the TLS energy splitting the reservoirs can excite the 
system and transport is dominated by incoherent se- 
quential emission and absorbtion processes. Except for 
one numerical simulation, 9 all previous investigations of 
spin-boson heat transporter— have considered only this 




FIG. 1: (Top) Schematic picture of thermal transport in a 
two-reservoir spin-boson model. A two-level system mediates 
energy transfer between two harmonic oscillator baths held 
at different temperatures. Possible experimental realizations: 
(center) phonons in two bulk leads couple through a nonlinear 
molecular junction, (bottom) electromagnetic fluctuations of 
two linear circuits couple through a superconducting qubit. 



regime. But when the bath temperatures decrease be- 
low the two-level splitting it is no longer possible to ex- 
cite the mediating system and sequential processes are 
blocked. In this case one needs to take into account co- 
herent two-boson processes which correspond to virtual 
excitation of the two-level system. The qualitative pic- 
ture is the same as in electron transport in quantum dots 
under Coulomb blockade. At low enough bias voltage the 
charging energy blocks exponentially the sequential tun- 
neling processes where individual electrons hop in and 
out of the dot incoherently^ The dominant transport 
processes then consist of coherent two-electron cotunnel- 
ing processes where the population of the dot remains the 
same after the process^ In analogy to electron trans- 
port, we use here the word "cotunneling" for the virtual 
two-boson processes giving rise to the low temperature 
heat conductance. By calculating two-boson cotunneling 
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rates we obtain an analytic description for thermal trans- 
port below the sequential-cotunneling crossover. Our re- 
sults are valid for weak bath coupling or for an arbitrary 
coupling strength when the noncommuting part of the 
two-level Hamiltonian and coupling operators is small. 
We present a detailed analysis of the crossover region 
and discuss the applicability of the sequential tunneling 
approximation. By employing the scaling equations of 
the Ohmic spin-boson model we can take into account 
rcnormalization of the spin-boson parameters. This al- 
lows us to complement the perturbative results by scaling 
effects and extract an anomalous temperature scaling of 
cotunneling thermal conductance. 

The paper is organized in the following way. In Section 
|n]we introduce the studied model in detail. In Scction lTTTl 
we calculate the cotunneling transition rates of relevant 
two-boson processes. We present a general formula for 
heat current from which one can obtain an analytic ex- 
pressions for the low-temperature and high-temperature 
limits and study the crossover behavior. In Section IIVI 
we use the scaling equations of the Ohmic spin-boson 
model to calculate temperature scaling of the cotunnel- 
ing thermal conductance and in Section [V] we conclude 
our analysis. 



II. MODEL 

We model our system with the usual spin-boson Hamil- 
tonian, separating the environmental oscillators into left 
(L) and right (R) baths: 
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Ho 
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u)i{a\ai + \), (2) 
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The TLS is described by level detuning e and transition 
amplitude A, while the baths are characterized by their 
spectral functions xk(w) = 27r^ igif Cf6(uj — u>i), K = 
L,R. Equations ([I])-© constitute the most generic two- 
bath spin-boson system^ 

It is convenient to diagonalize also the TLS part of Ho 
with a rotation of the pseudospin by an angle 9 in the 
x-z plane, with tan# = A/e. In this new basis we have 



ieL,R 



Ht = (o z cos 9 — a x sin ( 
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(4) 
(5) 



where ujo = V e 2 + A 2 . It should be noted that the valid- 
ity of our perturbation expansion requires that the TLS 
level broadenings r and Ti, given by Eqs. and (|T7|) 
below, are small compared to the level separation U3q. 
Since r ,i contain the couplings Cj as well as sin#, our 
results are valid either for weak coupling or when the 
noncommuting part of the coupling is small (9 <C 1). 



III. FORMULA FOR THE CURRENT 

Our calculation of heat current between the baths is 
based on the generalized Fermi Golden Rule^ 5 Given two 
eigenstates \i) and |/) of H , the perturbation Ht in- 
duces transitions between them at a rate 



r^ f = 2n\(f\T\i)\ 5(Ei-E f ), 



(6) 



where Ei and Ef are energies of the two states, and the 
scattering operator T is 

T = Ht + HtGoHt + HtGoHtGoHt ~\~ ■ ■ ■ (7) 
Here Go is the retarded propagator for H at energy E i: 

G ° = Ek-H + iv (8) 

We will consider only those processes where one quantum 
of energy is transported between the baths, i.e., the initial 
and final states arc related as 
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■■ a' K ,a K \i) 



(9) 
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where K and K' refer to arbitrary oscillators in different 
baths and iVj are their occupation numbers. In principle 
the initial and final states of the TLS could be different 
but below we show that these inelastic processes vanish 
for two-boson transport. Restriction to Eq. © forbids 
processes with an odd number of bosons and therefore wc 
only keep the even powers of Ht in Eq. ([7]), yielding 16 



T = H t GHt, 
where G is solved from a Dyson equation as 

1 



G = 



Ei — H — £ 



(10) 



(11) 



with self-energy S = HtGoHt- We consider the self- 
energy to the lowest order and keep only the imaginary 
part which regularizes the divergence at u>o- We will now 
show that in this case the off-diagonal elements of E can 
be neglected. Let us decompose the T matrix elements 
as 

(f\T\i)= Y, (mT\ka')(ka'\G\ja)(ja\H T \i), (12) 

where \j) and |fc) are states in the bath subspace, and \a) 
and \a') in the TLS subspace. Similarly we can write E 
as 



(ka' \HT\ma) (mer|-Hy| jo) 

Ei - E ma + IT) 

(k<7'\HT\mo)(mo\HT\jo) 



Ei — E„ 



IT) 



(13) 



where \o) is the state with the spin flipped with respect to 
|er). As noted above, we will only be concerned with the 
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imaginary part of the self-energy, and for the first term 
on the RHS it is nonzero when Ei = E mr7 . On the other 
hand, to lowest order the self-energy is used to regularize 
the divergence of sequential energy transfer, and in that 
case the coupling terms Ht in Eq. (fl"2j) produce energy- 
conserving processes, i.e., Ef — Eka> and Ei = Ej a . 
Thus we have E ma — Ej a , implying |m) = \j) and 
(ma\HT\jcr) — 0. The first term of Eq. (|13l) can therefore 
be neglected. For the second term a non-vanishing imagi- 
nary part requires Ei = E mS , Eq. ^ gives JSj = Ef, and 
the sequential condition is Et = E^i, which combine to 
E mS = Eka> ■ Therefore if a' = a we have |m) = \k) and 
the second term also vanishes. Only the spin-diagonal el- 
ements a' = a of the second term give a finite imaginary 
part. This also means that within our weak-coupling ap- 
proximation only spin-flipping processes, proportional to 
a x in Ht, contribute to S. 

Next we argue that the self-energy elements which are 
off-diagonal in the bath subspace can also be neglected. 
From Eq. ([13")) one can see that nonvanishing off-diagonal 



elements must have \k) = \j) , where a- 1 ' denotes 

either a creation or destruction operator for the baths. 
In this case the sum contains only two terms, with \m) = 

\j) or l m ) = a 2 li) • I n contrast, for diagonal elements 
\j) = \k) the sum goes over all states of the form |m) = 

o!m\j)- Since each term in the sum is proportional to the 
coupling, at least in our weak-coupling approximation the 
diagonal elements dominate and the off-diagonal ones can 
be discarded. Thus we obtain the expression 

-Im(jaj£|jcr) = \(ma\H T \ja)\ 6(E j(T - E mS ) 



,(t) 



gives 

To = T L e- UJo/TL +T R e~ UJa/TR , (16) 

Ti = T L +T R , (17) 

T K = ^n 2 e X K{^)[n K {^) + l}, K = L,R. (18) 

Here is the decay rate of the excited TLS into bath 
K with temperature Tk and corresponding Bose distri- 
bution UK- 

To derive a formula for the heat current we first write 
down an expression for heat absorbed from the left bath 
and emitted to the right bath, with er being the initial 
(and final) state of the TLS 



J 



(19) 



where the sum is over all possible initial states of the 
reservoirs (with a thermal distribution Wi), and all final 
states of the form © with K — fa an arbitrary oscillator 
in the left bath and K 1 = f R in the right bath. Energy 
conservation requires that ujf L = ojf R . Current in the 

other direction, J R \ rL , is defined analogously, and the net 

heat flow is J( CT ) = J { l\r-Jr\l which 

can be evaluated 

as 



sin 4 9 
2tt 



du LJXL{u)xR(u)[n L (uj) - n R (uj)] 



1 



1 



where upper signs refer to a = 0, lower signs to a = 1 



(20) 



(14) 



that is, the imaginary part of the self-energy element is 
half of the Golden- Rule transition rate out of state \ja) . 
Collecting everything yields 



2tt 



(f\H T \ja)(j*\H T \i) 



0,a 



Et — Ei, 



i 2 Fj a 



5{E f -Ei). (15) 



Evaluating the sum over the intermediate states \ja) is 
straightforward. We note that the calculation shows that 
even for cotunneling (with Ei ^ Ej a ) it is again only 
the spin-flipping terms, proportional to <r x in Ht, which 
give a non-vanishing contribution. The coupling term 
proportional to a z is therefore completely irrelevant for 
the two-boson weak-coupling current. This is intuitively 
plausible since such a term only produces a small shift in 
the system energy levels. 

Next we replace the lifetime TJ^ in Eq. (|15[) with the 
lifetime of a state where the TLS subspace is in state \a) 
but the bath subspaces are in thermal states. Then Tj a 
becomes independent of the bath state \ j) , and we call the 
two possible values r and r l5 corresponding to ground 
state and excited TLS, respectively. A short calculation 





FIG. 2: (Top) At temperatures comparable to the two-level 
energy splitting transport is dominated by incoherent sequen- 
tial processes. (Bottom) At low temperatures it is no longer 
possible to excite the two-level system and transport is dom- 
inated by coherent two-boson processes creating virtual exci- 
tations. 

We still have to combine the two partial currents 
to obtain the full heat flow through the system. First 
note that under steady-state conditions the probability 
for the TLS to be in ground state is Po = Fi/(r + Ti), 
and for the excited state Pi = 1— Po- Let us then split the 
currents into sequential and cotunneling contributions, 17 

= J S eq + J^ot i see a l so Fig. [21 For sequential pro- 
cesses the intermediate state has the TLS spin flipped, 
for cotunneling processes this flipping is only virtual and 
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does not contribute to Pq/\. Now, for sequential trans- 
port, one current-carrying process consists of two inco- 
herent tunneling events, and if we start with, say, the 
TLS in ground state, then the intermediate state is ex- 
cited. But this intermediate state can equally well be 
seen as the initial state of another current-carrying pro- 
cess. Because of this overlap, one must be careful to avoid 
double counting, and as further elucidated in Fig. [3J the 
full sequential current is J scq = P Q J^q = Pi Jscq- 



time 



r l"seq 



0"seq 



FIG. 3: Evolution of the TLS state due to sequential tunneling 
events. The current carried by these tunneling events can be 
calculated in two different ways. On one hand, we can start 
with the TLS in ground state, and then one current-carrying 
process flips the state twice, and the next process again starts 
from the ground state. The average current transported by 



these processes is J s( 



p r< > 

J "Scq ■ 



On the other hand, we can 



describe the same sequence of events by starting with the 
TLS in the excited state, giving J scq = PiJseq- Thus the 
two expressions for sequential current are not additive but 
alternatives to one another. 



On the other hand, cotunneling events are non- 
overlapping and therefore they contribute additively, 

Jcot — PqJcoI + PiJclt- Total heat current through the 
system, J = J scq + J cot , is therefore 



J = P J^ + P,J^ - J si 



(21) 



The sequential current can be calculated by taking the 
difference of the forward and backward Golden Rule rates 
at either junction, e.g., for the left junction, 



J scq = u)a(P T L e 



-ojo/T l 



P%T L ). 



(22) 



We can still simplify Eq. (|2T|) . As will be shown below, 
at high temperatures sequential tunneling dominates and 
all three terms on the RHS of Eq. d2~T|) arc equal in mag- 
nitude and thus the last two terms cancel. They can 
also be neglected at low temperatures since then Pi and 
J scq are exponentially suppressed. For these limits we 
therefore have 



J = F n J (0) . 



(23) 



Numerical calculations suggest that this expression is ac- 
tually valid within a few precent at all temperatures. 
Equation of the same form has been previously used by 
Flensbergi£ in the context of electron cotunneling. 



A. High-temperature limit 

Next we will derive the high and low-temperature lim- 
its for Eq. (|23[) . corresponding to sequential and cotunnel- 
ing, respectively. At high temperatures resonant energy 
transfer at lo = ujq dominates, and therefore the second 
energy denominator in Eq. (|20l) can be dropped. Con- 
tribution from the resonance peak can be calculated by 
taking the limit ^o,i 
^-5(u! — ljq). The current is then 



0, implying \u - 0J ± §r| 



J 



sin 2 uj XL(^o)XB,(^o)[nL(^o) ~ n R (u} )] 
Xl(w ) [2n L (w ) + l] + Xr(uq) [2ti r (uj ) + ] 



(24) 

This expression is equivalent to the sequential tunneling 
current calculated from Eq. (|22l) . and to the result de- 
rived by Segal and Nitzan.— 

The processes leading to Eq. correspond to se- 
quential absorption-emission of bosons as depicted in 
Fig. [21 In this limit the heat current is sensitive only to 
the value of the bath spectral functions at resonance. At 
temperatures Tl,r *C ujq expression AM]) becomes expo- 
nentially small, signalling the breakdown of the sequen- 
tial tunneling approximation. 



B. Low-temperature limit 

In the low-temperature limit, Tl r <C we can ignore 
the w-dependence of the energy denominators in Eq. (|20l) . 
In addition, the excitation of the TLS is exponentially 
suppressed, so that J = . Then the general expres- 
sion reduces to 



J 



2 sin 4 9 



efcjwxL(w)XiiM[™LM - n R (u})] .(25) 



Transport is dominated by two-boson cotunneling pro- 
cesses during which one boson is absorbed and one emit- 
ted as in Fig. [2] 

Let us remark that both Eqs. ([M]) and ([23)) can be 
cast in Landauer form, J — J du> ujT{uj) \til(uj) — n R (ia)\ , 
where T(w) is the transmission function. An impor- 
tant difference between the two cases is that in the se- 
quential regime T depends on the bath temperatures 
while for cotunneling the transmission is temperature- 
independent. One consequence is that the cotunneling 
current is symmetric with respect to the reversal of the 
temperature gradient, and therefore in contrast to the 
high-temperature case? 2 - the spin-boson system cannot 
act as a thermal rectifier at low temperatures. 

For concreteness we calculate the cotunneling current 
in the case where the bath spectral densities have a 
power-law form xk{^>) = Akui Sk , with < sk < 1 cor- 
responding to sub-Ohmic, sk — 1 to Ohmic, and sk > 1 
to super-Ohmic environment. Then we have 



J = 



2A L A R 



ITLUn 



sin^rMCM^-rg), (26) 
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where v = 2 + s L + sr, and T and £ are the gamma and 
zeta function, respectively. It is interesting to note that 
the temperature dependence of the current is sensitive 
only to the sum of the exponents Sl + sr and a precise 
nature of the individual baths is irrelevant. 



C. Crossover 

In Fig. H] we plot the numerical solution of Eq. (|23|) 
together with the limiting forms given by Eqs. (I24[) and 
(121)1) . We see that the crossover between sequential tun- 
neling and cotunneling is quite sharp, and the two ana- 
lytical formulas are very well applicable outside a narrow 
transition region. By patching together the two limiting 
expressions one achieves an accurate analytical descrip- 
tion of the transport problem at all temperatures. The 
crossover temperature depends on the system parameters 
logarithmically and it is generally of the order of 0.1 wq, 
below which the naive Golden-Rule treatment becomes 
invalid. 
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FIG. 4: Current through the system as a function of left- 
bath temperature Tl- Solid line is the numerical solution of 
Eq. (123[) , dashed line is the high-temperature approximation 
from Eq. (124[) , and dash-dotted line is the low-temperature 
approximation from Eq. (|26p. Both baths are ohmic, sl — 
sr = 1, with equal coupling strengths Al = Ar — 0.1, and 
sin0 = 1. Right-bath temperature is Tr = 0.95 Tl- 



IV. LOW TEMPERATURE PROPERTIES 
FROM RENORMALIZATION GROUP 



In the linear-response limit, T R -> T L = T, Eq. ([26]) 
shows that in the cotunneling regime thermal conduc- 
tance G depends on the temperature as G cx T 1+SL+SR . 
This result was derived using perturbation theory and 
as such relies on an assumption of weak coupling be- 
tween the TLS an the baths. In principle it is possible to 
calculate higher orders of perturbation theory and find 
corrections but this becomes soon unfeasible. Instead of 



carrying out a more complete perturbative calculation, 
in this section we use a renormalization group (RG) ar- 
guments to improve our analysis and find correction to 
the temperature exponent of G(T). 

To set up the RG framework, we introduce a cutoff 
frequency uj c above which the bath response functions 
Xl.r vanish. The initial value of lo c is irrelevant for our 
purposes. We then perform the RG step by lowering 
the cutoff and integrating out all bath oscillators with 
frequencies above the new cutoff. This system is again 
described by a Hamiltonian of the form ([TJ but with a 
new set of parameters. The process is continued until 
uj c is equal to the temperature T. This way the system 
parameters, which we have so far treated as constants, 
become effectively temperature dependent and give an 
extra contribution to G(T). 



0.5 




FIG. 5: Renormalization group flow of the parameters a and 
A/cue in the Ohmic spin-boson model, as given by Eqs. (|27[) 
and (|28|) . Two separatrices (thick lines) divide the flows into 
three qualitatively different regimes. 

The actual RG calculation is performed by the well- 
known Coulomb gas mapping A This procedure is only 
available for Ohmic baths, and therefore for the rest 
of this section we set s L = sr = 1. The partition 
function is transformed to a form which is identical to 
the anisotropic Kondo model after the Yuval- Anderson 
mapping*^ The renormalization procedure in this rep- 
resentation was worked out by Anderson, Yuval and 
Hamann in the seminal paper [191 ] . To conform with es- 
tablished notation we introduce dimensionless coupling 
constants (Xl,r such that the bath spectra are given by 
Xk{oj) = 7rcvifW. The RG flow equations for the one-bath 
spin-boson modeli^S are directly applicable for several 
baths, 



d(A/oj c ) 
d In lu c 

da 



-(1 



d In uj c 



UJ, 



(27) 
(28) 
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with a = ai + aji being a combined coupling for both 
baths. These equations are valid for a < 4 and A/uj c < 1. 
The flow of the parameters with decreasing lu c is depicted 
in Fig. [5j There a three qualitatively different regimes 
which we call A, B, and C. 

In regime A the renormalization of a is small in a real- 
istic case where the bare value of A/w c is small, so we can 
treat a as a constant. Then solving Eq. (f2"T)) yields^^ 

A = A (^V, (29) 
Wo/ 

where Ao and uj c q are the initial, unrenormalized values 
of A and ui c . The cutoff is taken all the way down to 
the temperature T, and we obtain A oc T a . Now we can 
again examine the T dependence of conductance. From 
Eq. (gfij), 

G oc Sin4 V - A4 T 3 oc [ T3+4Q ' 6 » A 
G a ~ (A 2 + e 2 )3 1 T 3 ~ 2a , e « A 

(30) 

Thus the RG calculation gives a non-perturbative correc- 
tion to the bare T 3 behavior depending on the interaction 
strength a. We have extracted two limiting forms of this 
correction, depending on the relative magnitudes of the 
TLS detuning and the renormalized tunneling element. 
Since A always renormalizes to zero, in the asymptotic 
low-temperature limit thermal conductance always obeys 
G oc T 3+Aa even if the bare parameters satisfy e< A. 

We would like to make a brief comment regarding the 
validity of the above considerations. There are two pos- 
sible restrictions, namely the validity of flow equations, 
Eqs. (|2T|) and (|28p. and the validity of the perturbation 
calculation for the current, Eq. (|2f)|) . Since a and A/oj c 
move towards smaller values, the flow equations remain 
valid throughout the process, all the way to oj c — > 0. 
Thus Eq. pop can be used even in the zero temperature 
limit. On the other hand, validity of Eq. (|26l) requires 
that the perturbation expansion parameter, essentially 
a sin 2 6?, is small. In regime A, where a > 1, this means 
that 8 must be small, i.e., A/e -C 1, so the possibility 
G oc T 3 ~ 2a is not realized in the region of validity of the 
perturbation theory and must be discarded. 

In regime B we can also accurately approximate a with 
a constant value, and therefore Eq. pU|) applies also to 
this case. However, now the flow takes A/uj c towards 
larger values eventually leaving the region of validity of 
the flow equations. Therefore the flow should be termi- 
nated when A/cj c = 1, and at this point the value of the 
cutoff isi20 

W c ,min — ^cO • 

Thus Eq. PH|) can only be used down to temperatures 
T ~ w c ,min- The perturbation parameter a sin 2 8 is al- 
ways small for any 8 by taking a small enough a, so 
both limiting cases of Eq. pU)) are relevant above the cut- 
off temperature w c . m i n . In the asymptotic weak-coupling 



limit a — > one recovers the simple G oc Independence 
as expected by the perturbation calculation. 

Regime C has the most complicated flow and a cannot 
be treated as a constant, prohibiting simple analytical 
estimates. Also in this case A/w c is a relevant perturba- 
tion eventually flowing to strong coupling. In principle 
the flow equations can still be used to analyze the tem- 
perature dependence of thermal conductance above the 
cutoff energy by numerically solving the flow. 

V. CONCLUSIONS 

We studied thermal transport in the spin-boson model 
at temperatures below the TLS energy splitting. At 
low temperatures a naive Golden- Rule approximation be- 
comes invalid and transport is dominated by two-boson 
processes. Qualitatively the behavior is similar to elec- 
tronic transport in quantum dots under Coulomb block- 
ade when lowest-order processes are exponentially sup- 
pressed. Applying perturbation methods we obtained an 
analytical description of low-temperature transport and 
investigated the crossover between the low-temperature 
and high-temperature limits. The crossover region is 
sharp and an accurate analytical description at all tem- 
peratures is provided by patching the two limiting cases 
together. In the case of reservoirs with simple power law 
spectral densities Xl/r( u ) ^ u Sl/r the low-temperature 
thermal conductance obeys G(T) oc T 1+SL+SR . Going 
beyond perturbation theory, we applied renormalization 
group arguments to extract anomalous temperature scal- 
ing of thermal conductance in the Ohmic case. The power 
law predicted by simple perturbation theory was modi- 
fied by interaction-dependent correction which becomes 
significant outside the weak-coupling regime. 

The studied problem is not only important for future 
applications, but also relevant in systems available for 
experimental studies presently. One of the most flex- 
ible systems to probe our predictions is a supercon- 
ducting flux qubit coupled to two linear circuits acting 
as reservoirs . 12 ' 21 In this realization the electromagnetic 
fluctuations of the reservoirs play the role of bosonic 
modes in the spin-boson model. Using external fields 
to tune the energy splitting of the qubit it is possible to 
address different regimes in a controlled manner. Con- 
sidering the generic role of the spin-boson model, it is 
expected that our results are important for a variety of 
applications in heattronics of molecular and nanoscale 
systems in the future. 

Upon completion of the manuscript we learned about a 
new work by Wu and SegaL 2 — They have studied a system 
similar to ours and their results agree with Eq. (|2"5|) . 
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